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ABSTRACT 

Recent theoretical investigations have pointed out that the cores of the hya lines of 
H i and He n should show measurable scattering polarization signals when observing 
the solar disk, and that the magnetic sensitivity, through the Hanle effect, of such lin- 
ear polarization signals is suitable for exploring the magnetism of the solar transition 
region. Such investigations were carried out in the limit of complete frequency redis- 
tribution (CRD) and neglecting quantum interference between the two upper /-levels 
of each line. Here we relax both approximations and show that the joint action of 
partial frequency redistribution (PRD) and /-state interference produces much more 
complex fractional linear polarization (Q/I) profiles, with large amplitudes in their 
wings. Such wing polarization signals turn out to be very sensitive to the temperature 
structure of the atmospheric model, so that they can be exploited for constraining the 
thermal properties of the solar chromosphere. Finally, we show that the approximation 
of CRD without /-state interference is however suitable for estimating the amplitude 
of the linear polarization signals in the core of the lines, where the Hanle effect oper- 
ates. 

Subject headings: polarization — scattering — radiative transfer — Sun: chromo- 
sphere — Sun: transition region — Sun: surface magnetism 



1. Introduction 

In the transition region between the chromosphere and corona of the Sun, the kinetic tempera- 
ture suddenly jumps from 10 4 K to 10 6 K and the plasma changes from partially to practically fully 



1 Institute) de Astroffsica de Canarias, E-38205 La Laguna, Tenerife, Spain 
2 Departamento de Astroffsica, Facultad de Ffsica, Universidad de La Laguna, Tenerife, Spain 
3 Consejo Superior de Investigaciones Cientfficas, Spain 
4 Astronomical Institute ASCR, Ondfejov, Czech Republic 



-2- 



ionized. Although various physical mechanisms have been proposed for explaining this tempera- 
ture increase, the lack of reliable measurements of key physical quantities, such as the magnetic 
field, represents today the most serious limitation to a better comprehension of the role played 
by this boundary region on the heating of the solar corona. In order to make measurements of the 
quantities that remain basically unknown, we need (1) to identify observables sensitive to the phys- 
ical conditions of the solar transition region, (2) to develop the instruments needed for measuring 
such observables, and (3) to infer the relevant physical quantities (e.g., the strength and orientation 
of the magnetic field) through realistic modeling of the measured observables. 



Recently, iTrujillo Bueno et al.1 (|201l|) argued that the hydrogen Lya line at 1216 A should 
show measurable scattering polarization signals when observing the solar disk, and that via the 
Hanle effect the line-center amplitudes of the Q/I and U /I linear polarization profiles must be sen- 
sitive to the strength and orientation of the magnetic field in the solar transition re gion (with good 



sensiti vity to magnetic field strengths between 10 and 100 G). In a subsequent paper. ITrujillo Bueno et al 



(120121) pointed out that significant line-center scattering polarization signals are to be expected also 
for the Lya line of He n at 304 A, and that for the magnetic field strengths expected at transition 
region heights (i.e., 5<100 G outside active regions) the He n 304 A line, due to its very large Ein- 
stein coefficient, is immune to the Hanle effect (i.e., it is a line whose observed linear polarization 
pattern could be used as a reference for facilitating the identification of the observational signature 
of the Hanle effect in the hydrogen Lya line). The possibility that scattering processes produce 
measurable linear polarization signals in these and other emission lines of the solar transition re- 
gion is of great scientific interest, because it is mainly through the Hanle effect that we may hope to 
explore the magnetism of the upper chromosphere and transition region of the Sun. As a matter of 
fact, with few exceptions, the polarization signals induced by the Zeeman effect in transition region 
lines are expected to be very weak (e.g., in typical semi-empirical models of the solar atmosphere 
a volume-filling magnetic field of 100 G inclined by 45° with respect to the line of sight (LOS) 
produces in both Lya lines circular polarization VI I signals significantly smaller than 0.1%, while 
the contribution of the transverse Zeeman effect to their linear polarization is insignificant). 



The conclusions of lTrujillo Bueno et al.l (1201 Ul2012|) were obtained through detailed radiative 

transfer (RT) calculations in semi-empirical and hydrodynamical models of the solar atmosphere, 

applyi ng the quantum theory of polarization described in the monograph by lLandi DeglTnnocenti & Landolfi 
(|2004|) . This density-matrix theory is based on the approximation of complete frequency redistri- 
bution (CRD), which neglects correlations bet ween the frequencies of the incom ing and outgoing 
photons in the scattering events. In addition, ITrujillo Bueno et al.l (|201 lL 120121) neglected quan- 



tum interference between the 2 Pi/ 2 and 2 P 3 / 2 upper levels of such Lya lines. Both approximations 
should be suitable for estimating the linear polarization amplitudes at the line center, which is 
where the Hanle effect operates. However, in strong resonance lines the joint action of partial 
frequency redistribution (PRD) and /-state interference can produce important observational sig- 
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natures in the wings of the fracti onal linear polarization (Q/I) profiles (e.g., Smitha et al. 20 11 al ibi: 
Belluzzi & Trujillo Buenoll2012l) . 



The aim of the present paper is to investigate the impact of PRD and /-state interference 
effects on the scattering polarization of the Lya lines of H i and He n. To this end, we apply 
the same theoretical framework and RT code that lBelluzzi & Trujillo Buenol (|2012l) developed for 
investigating the scattering polarization pattern across the h and k lines of Mg n. Like the Mg n 
h and k lines, also the above-mentioned Lya lines result from two transitions between a common 
lower level with angular momentum 7=1/2 and two upper levels with 7=1/2 and 7 = 3/2. 
However, in contrast with the case of the Mg n resonance lines, the two components of each Lya 
line are fully blended since the fine-structure splitting between the two upper /-levels is much 
smaller than the Doppler line widthfl 



2. Formulation of the problem 



This investigation has been carried out with i n the same theoretical framework that has been 
recently developed by iBelluzzi & Trujillo Buenol (|2012|) for investigating the polarization proper- 
ties of the Mg ii h and k lines. The reader is referred to that paper for more details on the theoretical 
tools that we have applied. 

We consider a two-term atomic model both for H i and He n, the lower term being composed 
of the ground level ( 2 Si/2) and the upper term of the upper levels ( 2 Pi/2 and 2 P3/2). Such atomic 
model allows us to take into account the fine structure of the Lya line, as well as the contribution 
of the /-state interference between the two upper /-levels. 

The Non-LTE RT problem is formulated within the redistribution matrix formalism. We con- 
sider a linear combination of the 7? n and R m redistribution matrices (which describe purely coherent 
scattering and completely redistributed scattering in the atom rest frame, respectively), for a two- 
term atom with unpolarized and infinitely sharp lower term, in the absence of magnetic fields. The 
redistribution matrix i? n in the ob server frame is calculated starting from the expression derived by 
Landi Pegl'Innocenti et al.1 (119971) in the atom rest fram e, and taking Doppler redistribution into 
account (see Eq. (3) of IBelluzzi & Trujillo Buenol 1201 21) . The redistribution matrix R m is calcu- 
lated in the limit of complete frequency redistribut ion (CRD) in the observer frame, s tarting from 
the CRD theory for a two-term atom presented in lLandi DeglTnnocenti & Landolfil (|2004|) . The 
effect of inelastic and superelastic collisions is included under the same approximations made in 



'Note that the fine-structure splitting between the two upper /-levels, though small, is in any case much larger than 
the levels natural width. 
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Belluzzi & Truiillo Bueno (2012). The collisional rates for the H i Lya line have b een calculated 
form lPrzybilla & Butler! (|2004|) . those for the He n Lyar line from ljanev et al.1 (|1987|) . 



We conside r the branching ratios for Ru and R m derived for the two-level atom case (see 



Bommier 1997a,b). Neglecting depolarizing collisions, such branching ratios are given by a and 



(1 - a), respectively, with@ 



r ff + r, 



a = 



r R + Y I + r E 



(i) 



The quantities Y R , Yj, and Y E are the broadening widths of the upper level due to radiative decays, 
superelastic collisions, and elastic collisions, respectively^ 

The variation of a as a function of height for the two li nes under investigation, in the semi- 
empirical solar atmosphere model C of Fontenla et al.1 (11993b (he reafter. FAL-C), is shown in the 
left panel of Fig. 1. As pointed out in lTrujillo Bueno et al.l (|2012|) . and as it can be seen from the 
right panel of Fig. 1, the He n Lyor line forms in a thin layer of the transition region (approximately, 
between about 2000 and 2200 km). As shown in the left panel of Fig. 1, at these heights a is 
essentially equal to unity so that the assumption of purely coherent scattering in the atom rest 
frame is an excellent approximation for the modeling of this line. The situation is slightly different 
as far as the H i Lya line is concerned. While the core of this line forms in the transition region, 
just a few kilometers below the core of the He n Lya line, its broad wings form deeper in the 
chromosphere, where a assumes values appreciably smaller than one, though larger than 0.9. The 
approximation of purely coherent scattering in the atom rest frame is thus not so good for this line, 
and its modeling requires to take both Ru and Rm into account. 

The polarization profiles of the two lines under investigation are calculated by applying a 
Non-LTE RT code based on the angle-averaged expressions of Ru and R m , and including the con- 
tribution of an unpolarized continuum. The numerical method of solution is a direct generaliza- 
tion t o the PRE) case of the Jacobian iterative scheme presented in lTrujillo Bueno & Manso Sainz 
(I1999|) . The initialization of the iterative calculation is done using the self-consistent solution of 
the corresponding unpolarized problem, obtained by applying Uitenbroek's (2001) RT code. The 
same code is used to calculate the continuum opacity (including the UV line haze) and emissivity. 

We point out that due to the negligible impact of /-state interference on Stokes-/, the PRD 
intensity 1(A) profiles obtained with our code coincide with those computed using Uitenbroek's 
(2001) code (which neglects /-state interference). As we shall see below, /-state interference 



2 Note that the branching ratios a and (1 - a), defined according to Eq. (1), do not include the factor (1 — e), with e 
the photon destruction probability. This factor is directly included in the redistribution matrices. 

3 Our theoretical formulation is based on the hypothesis that such broadening withs are identical for the two levels 
of the upper term. This is a very good approximation for the model atoms considered in this work. 
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efFects, together with those caused by PRD, have however an important impact on the Q/I profiles. 



3. The scattering polarization profile of the H i Lyar line at 1216 A 



Figure 2 shows the PRD Q/I profiles of the H i Lya line, obtained by taking into account 
and neglecting quantum interference between the two upper 7- levels. The calculations have been 
performed in the FAL-C atmospheric model, for a line-of-sight (LOS) with /i = cos 6 = 0.3 (with 
6 the heliocentric angle). 

The two profiles perfectly coincide i n the core of the line where 7-st ate interference does 
not produce any observable signature (cf., iBelfuzzi & Trujillo Buenoll201 In this spectral re- 



gion, our PRD profile show s a good agreement with the corresponding CRD profile calculated by 



Trujillo Bueno et al.1 (120111) (see the left panel of their Fig. 2) 



2SpC 

f 



Outside the line core, PRD effects produce a complex linear polarization profile with extended 
wings, while 7-state interference plays a significant role, producing much larger polarization am- 
plitudes with respect to the case in which it is neglected. The profile obtained by taking 7-state 
interference into account shows in particular two narrow negative peaks at approximately +0.4 A 
from line center, with an amplitude of about -7%, and two broad negative lobes with a minimum 
of about -6% at approximately +10 A from line center. As seen in Fig. 3, the shape of the Q/I 
profile remains qualitatively the same for line of sights corresponding to larger // values, although 
the Q/I amplitudes are clearly smaller. 

Figure 4 shows the sensitivity of the scattering polarization profile of the hydrogen Lya line to 
the thermal str ucture of the model atm osphere. We compare the results for the atmospheric models 
C, F and P of Fontenla et al.l (119931) . which can be considered as illustrative of quiet, network 
and plage regions. Such sensitivity is due to the fact that the anisotropy of the incident radiation 
field (which induces atomic level polarization in the 2 P 3 / 2 level and quantum interference between 
its sublevels and those of the 2 Pm level) depends on the gradient of the Stokes-7 comp onent of 
the source function (e.g., iTrujillo Buenol 1200 ll : LLandi DeglTnnocenti & Landolfil 120041) . Notice 
that especially the wing Q/I signals, which are insensitive to the magnetic field (since the Hanle 
effect operates only in the line core), are very sensitive to the temperature structure of the model 
atmosphere. Therefore, in addition to the 1(A) profile itself, the wings of the Q/I profile could help 
us to constrain the thermal properties of the observed atmospheric region, and this would in turn 
facilitate the modeling of the line-core signals and the inference of the magnetic field. 



4 We point out that our PRD calculations and the CRD calculations of lTrujillo Bueno et al.l (1201 11) were made fixing 
the total neutral hydrogen number density to the values tabulated in the FAL-C model. 
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4. The scattering polarization profile of the He n Lya line at 304 A 

The left panel of Fig. 5 shows our PRD results for the Q/I profile of the Lya line of He n 
calculated in the FAL-C atmospheric model, taking into account (solid line) and neglecting (dashed 
line) the impact of /-state interference. The line-core polarization, and its modification by the 
Hanle effect, can be safely modeled in the CRD limit, and neglecting /-state interference. However, 
as soon as we go outside the line-core region, the combined action of PRD and /-state interference 
effects produces strong Q/I signals, which are not obtained when such physical ingredients are 
neglected. The right panel of Fig. 3 shows the center to limb variation of the Q/I profile. 

Finally, in the right panel of Fig. 5 we provide information on the sensitivity of the Q/I profile 
to the thermal structure of the atmospheric model. As can be seen, the wings of the Q/I signals 
are very sensitive to the model's temperature structure, while the line-core signals are practically 
identical in the three atmospheric models. 



5. Understanding the impact of /-state interference 

The impact of /-sta te interference on the Q/I profiles of various multi plets was first demon- 
strated and explained in IStenflol (|1980|) and considered in greater detail in IStenflol (|1997|) . Landi 
DeglTnnocenti & Landolfi (2004), and belluzzi & Trujillo Buenol (boill) . 



The results presented in the previous sections show that in the wings of the lines, the amplitude 
of the Q/I profiles is about a factor 3 larger when the effects of /-state interference are taken into 
account (see Figs. 2 and 5). In order to understand this behavior, we first note that far away from 
its "center of gravity", a multiplet (in which the contribution of interference between different /- 
levels is taken into accou nt) behaves in resonance scattering as a two-level atom transition with 



J c = L { and /„ = L u (see lLandi DeglTnnocenti & Landolfi 



2004). We also observe that at these 



frequencies scattering is purely coherent in the observer frame. St arting from the expressions of 
the emission coefficients given in lLandi DeglTnnocenti et al.l (|1997|) . it can be shown that, far from 
the center of gravity of the multiplet, the fractional polarization of the radiation scattered at 90° is 
given by 



, ^ = -W 2 {L t ,L u )w{v) 

£/(v) 4 



(2) 



where the quantity W 2 (Jt, /«) is defined in Eq. (10.17) of lLandi DeglTnnocenti & Lan dolfi (2004), 
w{v) = V2/q(v)//q(v) is the monochroma tic anisotropy factor, and where we hav e neglected the 
contribution of /q to si (see Eq. (5.157) of lLandi DeglTnnocenti & Landolfi! 120041 for the defini- 
tion of the radiation field tensor /g(v)). 



From the same equations, it can be shown that for a multiplet having a single /-level in the 
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lower term, far from its center of gravity, the fractional polarization obtained taking into account 
the contribution of the various lines, but neglecting interference, is given by 

beWL int . = £ 7 ») w ^ ■ 0) 

where <S /„ is the relative strength of the J e — » J u transition (see, e.g., Eq. (3.65) of Landi Degl'Innocenti 
& Landolfi 2004). For our 2 5 - 2 P multiplet, having L e = and L u = 1, we then have 

beOOk ^2d,0) 



[Pfi(v)Ltot S l/2 W 2 (l/2,l/2) + S 3/2 W 2 (l/2,3/2) " 

Using the numerical values W 2 (1,0) = 1, W 2 (l/2,l/2) = 0, W 2 (l/2,3/2) = 1/2, S l/2 = 1/3, 
£3/2 = 2/3, we easily find [pg(v)] int / [Pe(v)] noint = 3, in agreement with the numerical results 
shown in this Letter. 



6. Conclusions 

As shown in this Letter, the joint action of PRD and /-state interference effects produce com- 
plex scattering polarization Q/I profiles in the Lyar lines of H i and He n, with large polarization 
amplitudes in the wings. PRD effects determine the qualitative shape of the Q/I profiles, while 
/-state interference plays a key role in producing the large polarization amplitudes in the wings of 
the lines. Such wing polarization signals turn out to be very sensitive to the temperature structure 
of the atmospheric model. They could thus be exploited, in addition to the intensity profile itself, 
to constrain the thermal properties of the solar chromosphere. 

As expected, the CRD approximation is only suitable for estimating the line-core signals. For 
instance, for a LOS with = 0.3 the line-center Q/I signal of the He n 304 A line calculated with 
PRD and 7-state interference in the FAL-C atmospheric model is about -1.3%, while it is about 
a factor three smaller for the hyd rogen Lyo- line. The s e line -center Q/I amplitudes are in good 



agreement with those obtained by iTrujillo Bueno et al.l (|2012f) . who assumed CRD and neglected 



7-state interference effects for investigating the impact of the Hanle effect on the linear polarization 
produced by scattering processes in the core of such Lya lines. 

Finally, it is of interest to estimate how large are the linear polarization amplitudes that result 
from the wavelength-integrated Stokes profiles ((Q) /(/)). It turns out that for the He n line (Q) /(I) 
is similar to the line-center amplitude (e.g., for a LOS with /i = 0.3 we have (Q)/(I)~ - 1.3% 
in the FAL-C model), while for the hydrogen Lya line (Q)/(I) is much larger than the line-center 
amplitude (e.g., for a LOS with /i = 0.3 we have {Q)/{I)~-3% in the FAL-C model). Narrow band 
filter polarimetry (e.g., with a FWHM~ 1 A) in these Lya lines would provide interesting /, Q/I 
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and U /I images over a large field of view. However, the information on the magnetic field of the 
solar transition region is encoded in the line-core polarization of the hydrogen Lya line, which is 
where the Hanle effect operates. For this reason, the Chromospheric Lya Spectropolarimeter (see 



Kobayashi et al.ll2012() aims at doing spectropolarimetry of the hydrogen Lya line with a spectral 



resolution of 0.1 A and a polarimetric sensitivity of 0.1%. Since outside active regions the He n 
304 A line (which originates entirely within the transition region) is practically immune to the 
Hanle effect, the possibility of having in addition the information provided by filter polarimetry 
in this reference line would facilitate the diagnostics of the thermal and magnetic structure of the 
solar transition region. 
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Fig. 1. — Left panel: the solid lines show, for each hya line, the variation with height in the 
FAL-C model atmosphere of the branching ratio a of Eq. (1). Note that for the He n 304 A line 
a is essentially equal to unity at chromospheric and transition region heights. The corresponding 
dotted lines give the atmospheric height where, at each wavelength, the ensuing optical depth is 
unity along a LOS with [i = 0.3 (see also the right panel for a magnified visualization). 
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Fig. 2. — The PRD Q/I profile across the hydrogen hya line, calculated in the FAL-C model 
atmosphere for a LOS with jj. = 0.3, taking into account (solid line) and neglecting (dashed line) 
/-state interference. The right panel shows in more detail the line core region. The reference 
direction for positive Q is the parallel to the nearest limb. 
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Fig. 3. — The center to limb variation of the Q/I profile of the Lya lines of H I (left panel) and 
He II (right panel), calculated in the FAL-C model atmosphere. The reference direction for positive 
Q is the parallel to the nearest limb. 




Fig. 4. — The Q/I profile of the hydrogen Lya line calculated in the indicated atmospheric models, 
taking into account PRD and 7-state interference effects. The right panel shows in more detail the 
line core region. The reference direction for positive Q is the parallel to the nearest limb. 
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Fig. 5. — Left panel: the PRD Q/I profile across the He n Lya line, calculated in the FAL-C model 
atmosphere for a LOS with n = 0.3, taking into account (solid line) and neglecting (dashed line) 
/-state interference. Right panel: the Q/I profile calculated in the indicated atmospheric models, 
taking into account PRD and /-state interference effects. The reference direction for positive Q is 
the parallel to the nearest limb. 



